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ABSTRACT 

The photometric and spectroscopic properties of 26 well observed Type la Supernovae (SNe la) 
were analyzed with the aim to explore SN la diversity. The sample includes (Branch-)normal SNe 
as well as extreme events like SNe 1991T and 1991bg, while the truly peculiar SN la, SN 2000cx 
(Li ct al. 2001) and SN 2002cx (Li et al. 2003) are not included in our sample. A statistical 
treatment reveals the existence of three different groups. The first group (FAINT) consists of 
faint SNe la similar to SN 1991bg, with low expansion velocities and rapid evolution of Sill 
velocity. A second group consists of "normal" SNe la, also with high temporal velocity gradient 
(HVG), but with brighter mean absolute magnitude < Mb >= —19.3 and higher expansion 
velocities than the FAINT SNe. The third group includes both "normal" and SN 1991T-like 
SNe la: these SNe populate a narrow strip in the Sill velocity evolution plot, with a low velocity 
gradient (LVG), but have absolute magnitudes similar to HVGs. While the FAINT and HVG 
SNe la together seem to define a relation between 7^(SiII) and Arnis(B), the LVG ones either do 
not conform with that relation or define a new, looser one. The 7?.(SiII) pre-maximum evolution of 
HVGs is strikingly different from that of LVGs. The impact of this evidence on the understanding 
of SN la diversity, in terms of explosion mechanisms, degree of ejecta mixing, and ejecta-CSM 
interaction, is discussed. 



Subject headings: supernovae: general 

1. Introduction 

Given the role of Type la Supernovae as dis- 
tance indicators for cosmology and main producers 
of heavy elements in the Universe, understanding 
the physics of their explosions and how it influ- 
ences observables is one of the most fundamental 
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issues in modern Astrophysics. One of the keys to 
penetrate into the secrets of SN la physics is to 
explore the origin of their diversity. 

During the last decade a new paradigm for 
SNe la was developed. In particular, a correla- 
tion between the peak luminosity and the shape 
of the early light curve was found, with brighter 
objects having a slower rate of decline than dim- 
mer ones (Phillips 1993; Phillips et al. 1999, 
hereinafter P99). This is matched by a spec- 
troscopic sequence, defined by the ratio of the 
depth of two absorption features of Sill at 5972 
and 6355 A (typically observed at 5800 and 6150 
A, respectively) (Nugent et al. 1995). This ratio, 
7^(SiII), also correlates with the absolute mag- 
nitude of SNe la and, in turn, with the rate of 
decline. Spectral modeling indicates that most of 
the spectral differences are caused by variations in 
the effective temperatures which, in the context 
of Chandrasekhar-mass explosions, can be inter- 
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preted in terms of a variation in the mass of 56 Ni 
produced in the explosions. 

Alternatively, Garnavich et al. (2004), us- 
ing synthetic spectra, tentatively suggest that 
in SN 1999by and, in general, in SNe la with 
Ami5(B)> 1.2 the 5800 A feature is mostly due 
to Till transitions rather than to Sill. 

Although a one-parameter description of SNe 
la has proved to be very useful, it does not 
completely account for the observed diversity of 
SNe la, (e.g. (Benetti et al. 2004; Pignata et al. 
2004)). Hatano et al. (2000) showed that, contrary 
to expectations, 7^(SiII) correlates poorly with 
the photospheric velocity deduced from the Si II 
A6355 absorption. To account for this, they sug- 
gested that two or more explosion mechanisms are 
required to explain SN la variety. Furthermore, 
some SNe la with normal spectra (i.e. showing 
only lines of typical ions) were noticed to show ex- 
ceptionally high absorption line blueshifts (Branch 
1987). Finally, no correlation was found between 
the blueshift of the Si II 6355 absorption at the 
time of maximum brightness and the decline rate 
parameter, Ami5(B), in a small sample of well- 
observed SNe la (see also Patat et al. 1996). 

In this paper we further explore this issue, using 
detailed observations of a large sample of SNe la. 

2. Analysis 

The sample used in this work consists of 26 well- 
studied SNe la. It includes (Branch-)normal SNe 
(Branch & van den Bergh 1993) as well as extreme 
events like SNe 1991T and 1991bg. Photometric 
parameters such as Arai5(B), extinction and ap- 
parent magnitude at maximum are from P99 or al- 
ternatively from Altavilla et al. (2004) . The spec- 
tral parameters have been (re-)measured homoge- 
neously on the available spectra, most of which 
are published although unpublished material col- 
lected in the Asiago Supernova Archive (ASA) by 
the Padua- Asiago SN group and the European Su- 
pernova Collaboration was also used. In Table 1 
we summarize the main photometric and spectro- 
scopic parameters of the SN la sample, together 
with the morphological type (T) of the host galaxy 
as given in RC3. The SNe are divided into three 
groups according to the criteria discussed in Sect. 
2.1. In each group the SNe are arranged by de- 
creasing v (see Sect. 2.1). 



2.1. Expansion velocities from Sill 

The expansion velocity of the ejecta gives a di- 
rect indication of the kinetic energy of the explo- 
sion. The blue-shift of the Sill A6355 absorption, 
the most prominent line in the photospheric phase 
spectrum, traces the evolution of the expansion 
velocity of the ejecta over the first 4-5 weeks past 
explosion. A plot of the time evolution of the ex- 
pansion velocity, v(Si) for the SNe la of our sample 
(Figure 1) confirms the results of Branch (1987). 
At any given phase, expansion velocities span a 
wide range of <~ 4000 km s _1 , from the very rapidly 
expanding ejecta of SNe 1983G and 1997bp to the 
low velocity, SN 1991bg-like events. Between these 
limits both the "Branch normal" SNe la and the 
luminous SN 1991T-like events are found. A con- 
venient parameter to distinguish different SNe la 
is the expansion velocity measured 10 days past 
maximum, wi (SiII). As we mentioned earlier, this 
parameter correlates poorly with 72.(SiII) (Hatano 
et al. 2000) or with Am 15 (B) (Patat et al. 1996), 
hence it may give interesting clues for origin of 
SN la diversity. 

Upon a more careful scrutiny, and thanks to 
the excellent coverage offered by the new data, it 
can be noticed that the velocities show not only 
a spread in value, but also significant differences 
in velocity evolution. In particular, a sub-group 
of very homogeneous SNe la, characterized by a 
shallow evolution of the expansion velocity (open 
symbols in Fig. 1), populates a narrow strip in the 
post-maximum velocity evolution diagram. For a 
quantitative analysis we introduce a new parame- 
ter, v = — Av/At, which is the average daily rate 
of decrease of the expansion velocity (Tab. 1, col. 
4). This is derived from least squares fits of the 
measurements taken between maximum and either 
the time the Sill feature disappears or the last 
available spectrum, whichever is earlier. 

On average, the low (temporal) velocity gradi- 
ent (LVG) SNe la have a velocity gradient v < 
60 - 70 kms" 1 d" 1 . The SNe la shown as filled 
symbols have a larger v > 70 kms _1 d _1 , reach- 
ing 110-125 kms-M- 1 for SNe 2002bo, 1991bg, 
and 1983G (although the value for SN 1983G is 
quite uncertain). 

A similar dichotomy in velocity slopes between 
SNe la can be found in the velocities deduced 
from the SII A5640 line, this time at pre-maximum 
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Fig. 1. — Time evolution of the photospheric velocity derived from Sill A6355 for the SNe of Tab. 1. 
filled symbols refer to HVG, open symbols to LVG and starred to FAINT SN la. 
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phases (see Figure 11 of Benetti et al. (2004)). 

For an objective identification of possibly ho- 
mogeneous groups, we performed a hierarchical 
cluster analysis for the SN la sample of Table 1, 
considering both photometric and spectroscopic 
parameters, but neglecting the dependence on 
galaxy morphological type. Hierarchical Clus- 
ter Analysis (Anderberg 1973) is an exploratory 
data analysis tool which aims to identify relatively 
homogeneous groups of events based on selected 
characteristics, using an algorithm that starts with 
each case in a separate cluster and combines clus- 
ters until only one is left. Cluster analysis sim- 
ply discovers structures in data without explain- 
ing why they exist. The choice of the number of 
clusters to be considered is somewhat arbitrary al- 
though a criterion is the distance of the groups in 
the parameter space. 

Indeed, in the five dimensions space we found three 
well separated clusters: 

a) A first cluster (FAINT) consists of faint 
SNe la, with < M B >= -17.2 (a = 0.6), 
similar to SN 1991bg. All these SNe 
have a high post-maximum luminosity de- 



cline rate, <Ami 5 (B)>= 1.83 (a = 0.09), 
and Sill line ratio, <7e(SiII)>=0.58 (cr = 
0.05). They have small expansion veloci- 
ties, <«io(SiII)>=9.2 (cr = 0.6), and a large 
velocity gradient, < v >= 87 (cr = 20). 

b) A second group consists of "normal" SNe la 

with high velocity gradient (HVG), < v >— 
97 (cr = 16). These SNe have average abso- 
lute magnitude < Mb > = —19.3, (cr = 0.3), 
Am 15 (B)=1.2 (cr = 0.1), and <^(SiII)>=0.20 
(cr = 0.05). They typically have high expan- 
sion velocities, <Vi (SiII)> = 12.2 (cr = 1.1). 

c) A third group consists of SNe with a low veloc- 

ity gradient (LVG): < v >= 37 (cr = 18), but 
it includes both "normal" SNe la and all the 
brightest SNe. Although its post-maximum 
decline rate is somewhat slower than that 
of the HVG SNe, <Ami 5 (B)>=l.l (cr = 
0.2), its average luminosity is similar: < 
M B > = -19.2(cr = 0.3), and so is their 
<7£(SiII)>=0.25 (cr = 0.07). On average, 
the LVG SNe have lower and more homo- 
geneous expansion velocities than the HVG 
SNe, <uio(SiII)> = 10.3 (cr = 0.3). 
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Fig. 2. — The pre-maximum temporal evolution of 
the 7?.(SiII) parameter for our sample of SNe la. 
Symbols are as in Fig. 1. 

For a few events, in particular SNe 1989B and 
1992A, cluster membership is uncertain. A small 
variation of the parameters, still within the errors, 
could shift them from the LVG to the HVG group. 

It is also interesting that on average the host 
galaxy morphological type is progressively later as 
one moves from FAINT to HVG and to LVG. The 
average values are —1.4, 0.6, and 2.5, respectively. 
The dispersion within each group is however very 
large (cr = 3 — 4). This is consistent with the find- 
ing that bright SNe la occur preferentially in late 
type galaxies, while faint SNe la are more often 
found in early type galaxies (Altavilla et al. 2004). 

2.2. Pre-maximum evolution of TvL(Sill) 

Figure 2 shows the pre-maximum evolution of 
7?.(SiII) for our sample of SNe la. Interestingly, 
the HVG SNe la for which very early observations 
are available show a dramatic temporal evolution 
of 7?.(SiII), starting from a high value well before 
maximum and leveling-out just before maximum, 
as was the case for SN 2002bo (Benetti ct al. 2004). 
On the other hand, the LVG SNe show on aver- 
age either no evolution in 72.(SiII) before maximum 
or an evolution in the opposite sense in the case 
of SN 1990N. Clearly, the number of SNe la for 
which very early spectra are available is still too 
small to draw definite conclusions, but it is worth 



speculating on a possible interpretation. 

Since 7£(SiII) is related to the photospheric 
temperature of the ejecta (Nugent et al. 1995), the 
curves in Fig. 2 should trace the temperature evo- 
lution in the line forming regions before maximum 
light. HVG SNe seem to start at cooler tempera- 
tures, which then increase approaching maximum. 
LVG SNe, on the other hand, have high temper- 
atures already well before maximum. This may 
be related to the fact that LVG SNe have lower 
expansion velocities, especially before maximum, 
and are also hotter. 

2.3. n(SHI) max vs. Am 15 (B) 

In Figure 3a, the value of TZ(Sill) ma x mea- 
sured for each SN la at maximum light is plot- 
ted against Amis(B). SNe in both the FAINT 
and HVG groups seem to follow the relation be- 
tween lZ(SilT) max and Amis(B) established by 
Nugent et al. (1995). LVG SNe, on the other hand, 
are more scattered in this plot, especially at the 
bright, slow end. Either they do not conform with 
the TZ(SilT) max - Ami5(B) relation or they define 
a new, looser one. 

The scatter of LVG SNe in Figure 3a, especially 
at low Ami5(B), suggests that another physical 
parameter besides the temperature is needed to 
describe their behavior. Garnavich et al. (2004), 
using SYNOW synthetic spectra, tentatively sug- 
gest that in SN 1999by and, in general, in SNe la 
with Ami 5 (B)> 1.2 the 5800 A feature is mostly 
due to Till transitions rather than to Sill, and 
claim that only with this interpretation can the 
lZ(SilT) max - Ami5(B) relation be understood. 
However, neither our synthetic spectral analysis 
of SN 1991bg (Mazzali et al. 1997), nor the spec- 
tral tomography analysis of the normal SN la, 
SN 2002bo exploring a wide range of photospheric 
temperatures and chemical compositions (Stehle 
et al. 2004), requires a relevant contribution from 
the Till transitions to fit the 5800 A feature. In- 
deed, in the above mentioned models, the 5800 A 
feature is always well fitted by the Sill 5972 A 
transition. 

As for the SNe la with Ami5(B)< 1.2, Garnavich 
et al. (2004) confirm that the 5800 A feature is 
indeed Sill. 
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Table 1: Observed parameters of SN la sample, the host galaxies and references. 



SN 


Arais(B)* 




V 


fio(Sill)*" 
[kms _1 1000 _1 ] 


K(SiII)™ x * 


T(KC3) 


References 








[kms _1 d _1 ] 














LVG 










92A 


1.47 ± 0.05 


-18.81 


67 ± 7 


10.83 ± 0.15 


0.38 ± 0.05 


-1.9 


P99; ASA; K93 


89B 


1.34 ± 0.07 


-18.87 


66 ± 5 


10.03 ± 0.10 


0.29 ± 0.05 


3.0 


P99; B90; W94 


03kf 


1.01 ± 0.05 


-19.37 


50 ± 5 


10.63 ± 0.10 


0.18 ± 0.05 


3.0 


a 


96X 


1.25 ± 0.05 


-19.24 


46 ± 5 


10.68 ±0.15 


0.25 ± 0.05 


-5.0 


P99; S01 


99ee 


0.94 ± 0.04 


-19.46 


42 ± 5 


9.90 ± 0.10 


0.22 ± 0.05 


4.0 


S02; H02 


90N 


1.08 ± 0.05 


-19.23 


41 ± 5 


10.12 ± 0.10 


0.21 ± 0.05 


3.8 


P99; ASA; L91 


94D 


1.32 ± 0.05 


-19.06 


39 ± 5 


9.89 ±0.10 


0.33 ± 0.05 


-2.0 


P99; P96 


03du 


1.06 ± 0.06 


-18.93 


31 ± 5 


10.10 ± 0.10 


0.22 ± 0.03 


8.0 


b 


Old 


1.15 ± 0.04 


-18.71 


31 ± 5 


10.24 ± 0.10 


0.30 ± 0.07 


5.9 


K03; W03; c 


97br 


1.04 ± 0.15 


-19.62 


25 ± 5 


10.40 ± 0.10 


< 0.1 f 


7.0 


ASA; L99 


99cw 


0.94 ± 0.05 


-19.24 


22 ± 7 


10.79 ± 0.15 


< 0.31 n 


1.5 


d 


91T 


0.95 ± 0.05 


-19.62 


11 ± 5 


9.87 ±0.10 


0.14 ± 0.05 


3.8 


P99; P92; M95 


98bu 


1.04 ± 0.05 


-19.12 


10 ± 5 


10.29 ± 0.10 


0.21 ± 0.05 


2.0 


P99; ASA; H00 








HVG 










83G 


1.37 ± 0.10 


-18.62 


125 ± 20 


13.49 ± 0.20 


0.28 ± 0.06 


-2.2 


H83; B85; T85; B91; B89; e 


02bo 


1.17 ± 0.05 


-19.42 


110 ± 7 


11.73 ± 0.15 


0.17 ± 0.05 


1.0 


B04 


97bp 


1.09 ± 0.10 


-19.69 


106 ± 7 


13.87 ±0.20 


0.16 ± 0.05 


# 


A04; ASA 


02cr 


1.33 ± 0.04 


-19.45 


92 ± 5 


10.52 ± 0.10 


0.23 ± 0.05 


1.0 


P04; f 


84A 


1.21 ± 0.10 


-19.46 


92 ± 10 


12.64 ± 0.15 


0.23 ± 0.05 


1.0 


B89; Ba89; W87; e 


89A 


1.06 ± 0.10 


-19.21 


90 ± 10 


12.13 ± 0.15 


0.22 ± 0.07 


4.1 


B91; e 


02dj 


1.12 ± 0.05 


-19.05 


86 ± 6 


12.18 ± 0.10 


0.17 ± 0.05 


-5.0 


g 


81B 


1.11 ± 0.07 


-19.21 


76 ± 7 


11.11 ± 0.15 


0.16 ± 0.05 


4.5 


P99; B83 








FAINT 










99by 


1.87 ± 0.10 


-16.64 


110 ± 10 


9.06 ± 0.10 


0.58 ± 0.06 


3.0 


B99; V02; H01 


91bg 


1.93 ± 0.10 


-16.81 


104 ± 7 


8.50 ± 0.15 


0.62 ± 0.05 


-4.7 


P99; T96 


97cn 


1.86 ± 0.10 


-16.95 


83 ± 10 


8.81 ± 0.20 


0.63 ± 0.06 


-5.0 


T98 


93H 


1.70 ± 0.10 


-18.20 


73 ± 8 


10.10 ± 0.20 


0.52 ± 0.05 


1.9 


P99; ASA; h 


86G 


1.78 ± 0.07 


-17.48 


64 ± 5 


9.36 ± 0.15 


0.53 ± 0.05 


-2.2 


P99; P87; C92 



* reddening corrected according to P99 



** Cepheids distances taken from Altavilla et al. 2004 (Ho = 72kms~ 1 Mpc -1 ) when available, otherwise relative 
distances to Virgo taken from Kraan-Korteweg (1982) and a Virgo distance of 15.3 Mpc (Freedman et al. 2001) or 
from Hubble flow 
*** in units of 1000 kins" 1 
**** measured at maximum light 

fSill A5972 not visible on the -3 d spectrum; spectra close to maximum not available, 
ffmeasured on the earliest available spectrum(</> = +5d) 
# Irregular galaxy 

a) Salvo et al. 2004, in preparation 

b) Stanishev et al. 2004, in preparation 

c) Mattila et al. 2004, in preparation 

d) Bufano et al. 2004, in preparation 

e) McDonald archive; 89A, 84A and 83G: determination of Ami 5 (B) from published data 

f) Kotak et al. 2004, in preparation 

g) Pignata et al. 2004, in preparation 

h) CTIO Archive 
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3. DISCUSSION 




Am 15 (B) 

Fig. 3.— a). Amis(B) vs. TZ(Sill) max for the SNe 
of our sample, b). Amig(B) vs. v of Sill A6355. 
Open symbols refer to LVG SNe, filled symbols to 
HVG SNe, and starred symbols to FAINT SNe. 

2.4. v vs. Ami 5 (B) 

The Amig(B) parameter is plotted versus v 
in Figure 3b. The expansion velocity evolution 
gradient, v, seems to be weakly correlated with 
Amis(B): while SNe with a large Amis(B) (the 
FAINT group) have a large v, normal SNe can 
have both a large or a small v. Cluster anal- 
ysis however suggests that we are dealing with 
three distinct families of SNe la: FAINT, LVG, 
and HVG. These three groups may be character- 
ized by different physical parameters governing the 
same explosion mechanism (possibly distinguish- 
ing LVG and HVG) or by a totally different kind 
of explosion, which may be the case especially for 
the FAINT group. 



Based on the evidence presented above, we can 
make a preliminary attempt to explain the causes 
of the spectroscopic diversity among SNe la. 

As Figures 1,2 and 3 show, SNe la can be di- 
vided into three groups. Each group has distinct 
physical properties, different from those of other 
groups. In particular, the FAINT group (essen- 
tially the SN 1991bg-like objects, plus SNe 1986G 
and 1993H according to this method - but see 
Meikle (2000) for an infrared view) clearly differs 
from the other two: these SNe are fast decliners in 
both luminosity and velocity, they have typically 
low velocities and occur in earlier- type galaxies. 
This may be not surprising, since in many ways 
SN 1991bg-like objects stand out as odd. 

For the other two groups, the situation is more 
complicated: both HVG and LVG include normal 
SNe, but the LVG include also all the brightest, 
slowest declining SNe. Our analysis suggests that 
LVG and HVG are two distinct groups, but they 
may possibly represent a continuum of properties. 

Interestingly, although it is common to refer to 
two main groups of peculiar SNe la, SN 1991bg- 
and SN 1991T-like, the latter SNe fall in the same 
class as normal SNe, while the former do not. 
There appears to be a discontinuity of proper- 
ties between SN 1991bg-like objects and all other 
SNe la, which is not seen for SN 1991T-like SNe. 
This is not what we might expect if SNe la be- 
haved as a simple one-parameter family of events. 

Based on a qualitative analysis Benetti et al. 
(2004) suggested that the large blueshift of the Sill 
line in SN 2002bo (a HVG SN) maybe the result of 
a delayed detonation explosion. Moreover, Lentz 
et al. (2000) find, from detailed non-LTE calcu- 
lations, that some delayed-detonation models can 
account for the very high Sill blueshift of another 
HVG, SN 1984A. It may be that the HVG SNe 
are delayed-detonations, their internal dispersion 
arising from a range of transition densities, while 
the LVG SNe are deflagrations. This would be an 
extension of the results of Hatano et al. (2000) 
(see their Figure 1). The only difference with 
their conclusions would be that SNe 1981B and 
1992 A would become a delayed-detonation and a 
deflagration event, respectively. From spectropo- 
larimetric studies Wang et al. (2004) also reached 
the conclusion that the explosion mechanism of 
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SNe 1984A, 1997bp, 2002bo (all HVG SNe) and 
of SN 2004dt, which is most probably also a HVG 
SN, may be markedly different from that of lower 
velocity objects such as SN 1994D. 

The LVG include all three SNe in our sam- 
ple with Ami 5 (B)< 1: SNe 1991T, 1999cw and 
1999ee. SN 1991T has also often been discussed 
as the result of a delayed detonation, especially in 
order to explain the high abundance of 56 Ni and 
its decay products at the highest velocities (Maz- 
zali, Danziger, & Turatto 1995), so its inclusion in 
the LVG does not support the hypothesis that all 
LVG are deflagrations. Very early measurements 
of v(Si) are not possible for SN 1991T, since the 
Sill line was almost absent in the earliest spectra 
owing to the high degree of ionization. The pho- 
tospheric velocities inferred from spectral models 
were however very high (Mazzali, Danziger, & Tu- 
ratto 1995). This might support the suggestion 
made by Wang et al. (2004) that SN 1991T-likc 
events could be 1984A-like events viewed at dif- 
ferent angles (and thus HVGs). 

Differences in the properties of the outer ejecta, 
such as a different degrees of mixing, or of circum- 
stcllar interaction, may also be at the origin of the 
difference between LVG and HVG. 

More efficient mixing out of heavy elements 
might result in an initially higher photosphcric 
velocity. At the earliest times the photosphere 
should in fact tend to trace the heavier elements, 
since they have much larger line opacity and this is 
the major contributor to the optical depth (Paul- 
drach ct al. 1996). The initial rate of decrease 
of the velocity with time would consequently be 
larger, as the photosphere moves inward to layers 
that are not so different from those of less mixed 
SNe. This would also mean an initially lower tem- 
perature (resulting from the large photosphere) 
but increasing with time. This would be the HVG 
group. 

On the other hand, less efficient mixing could 
lead to initially smaller photospheric velocities: 
light elements contribute much less to the opac- 
ity and thus the photosphere would be deeper and 
the pre-maximum temperature higher. The de- 
cline rate of the velocity would then be smaller, 
and the temperature would either stay constant 
or decline, depending on the exact combination of 
increasing luminosity and decreasing photosphcric 
velocity. These are the properties of the LVG. 



Very high-velocity features have been observed 
in all SNe with sufficiently early spectra (Mazzali 
ct al., in preparation). Maybe in the LVG the 
interaction affects the spectra only very early, as 
in SN 1999ee (Mazzali et al. 2005) or SN 1990N 
(Fisher et al. 1997; Mazzali 2001), and there is a 
sudden drop to lower velocities when the interac- 
tion ends, resulting in a lower post-maximum v, 
while it continues for a longer time in the HVG, 
so that wio(Sill) spans a broader range of values. 

In both of these last scenarios, LVG and HVG 
SNe would not necessarily be differentiated by the 
nature of the explosion, and they may even repre- 
sent a continuum of properties. 

We are in debt to J. Sollerman, P. Lundqvist, 
M.M. Phillips, S. Mattila and J.C Wheeler for 
providing us unpublished SNe spectra. We also 
acknowledge D.A. Howell and L. Wang for sending 
us the spectra of SNe 1999by and 2001el respec- 
tively. This work is supported in part by the 
European Community's Human Potential Pro- 
gramme under contract HPRN-CT-2002-00303, 
"The Physics of Type la Supernovae" . 
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